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Abstract 

In view of the mechanics problems of the thick-walled open-sectional structural members that cannot be dealt with by existing 
beam theories, a novel computational methodology named Finite Nodal-Line Method (FNLM) is developed to perform the 
analysis of the thick-walled open-section beam with a length-aspect ratio less than 3. The computational results coming from an 
example of L-shaped thick-walled beam have verified the validity of the FNLM. 
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Introduction 

With the rapid development of Chinese cities, almost all civil buildings are now built to be tall. In order to expand 
the utilization of the construction area, and to improve the indoor visual effect and convenient decoration, some 
new structural systems are widely adopted by the tall buildings such as a frame system with special-shaped 
columns, a frame shear wall system with special-shaped columns, and etc.. The shape of the cross sections of the 
special-shaped columns is complicated and the size is not small, some of their limb length close to or more than 1 
meters, so that the slenderness ratio (shear span ratio) of the special-shaped column at around 3 or less than 3. 

Consequently, a formidable challenge comes out, i.e., the application conditions for all current beam theory do not 
include tht thick-walled open section beam with a length-aspect ratio less than 3. In addition, due to the 
inconsistency between the torsional center and geometric centroid of a complex cross section such as channel 
shaped column, L shaped column, etc., a torsional deformation will be produced under the action of a transverse 
force without passing through the torsional center. For reinforced concrete material components, the warping 
normal stress caused by restrained torsion must be particularly concerned. It is just the fatal weakness of existing 
beam theory, i.e., when a short and thick- walled structural member of complex cross section cannot be freely 
twisted, the warping stress on the cross section is too hard to be obtained, although there are many beam theories 
that concern the torsion of a member [1-5]. 

In order to cope with the challenge, a novel computational methodology named Finite Nodal-Line Method 
(FNLM) is developed to perform the analysis of the thick-walled beam with a length-aspect ratio less than 3. 

Finite Nodal-line Method (FNLM) 

The displacement field of a structural member is composed of two parts: the first part is the longitudinal (or axial) 
displacement, and the second part is the transverse displacement. All the deformed cross sections, which are 
curved surfaces, are defined to be the longitudinal displacements, and all the torsional rotation angles of each cross 
section as well as the point set located on the deformed axial line are defined as the transverse displacement of the 
member. The longitudinal displacement and the lateral displacement together constitute the displacement field of 
the member. 

Without loss of any generality, the displacement field of a thick- walled member can be described by formulating 
the displacement field of a thick- walled beam with channel cross section as shown in Fig. 1(a). 
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FIG.l: A CHANNEL-SECTION CANTILEVER BEAM 

(a) NODAL LINES 

(b) STRESS STATE AT AN ARBITRARY POINT 


The nodal lines are illustrated in Fig 1(a). If a finite number of unknown functions ( vg(x) , / = 1,2,- • -n ) representing 
the longitudinal displacements of "n", nodal lines are selected to be basic unknown functions, and interpolation 
functions (<p i (y i ,z i ),i = l,2,---,n ) are constructed in between the nodal lines, the longitudinal displacement of the 
beam can be expressed as follows. 

u(x,y,z) = 'Zw i (x)<P i (y,z) = <P r (. y,z)W(x ), (1) 


in which. 
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where W(x) is the unknown function array of longitudinal displacements, and @(y, z) is the interpolation function 
array. The domain of w i (x) is the set of points located along the x-axis; the domain of each <p i (y i ,z i ) is a part area 
of the cross section enclosed by the connection lines of the node i and its nearest neighbors, which are the 
intersection points between the cross section and nodal lines. 

Similarly, if three other basic unknown functions v 0y (x) , v 0z (x) , and 6{x) are chosen to represent the transverse 

displacements of the beam axis in the y-z plane and the rotation around the axis x, such functions can be expressed 
as 




V(x) = \ 


d(x) 


( 3 ) 


The domain of the three unknown functions is the point set located on the beam axis. 

By means of the displacement field defined in equations (1) and (3), the displacement of an arbitrary point in the 
beam in the x, y, and z directions can be expressed as follows. 


u(x, y,z) = <P T (y,z)W(x) 
v y (x,z) = v 0y (x)-z&(x) >. 

v. (x, y) = v 0 _ (x) + y0(x) 


( 4 ) 


As a result, the state of strain at an arbitrary point in the beam can be expressed in terms of equation (4) as 
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In the above, the comma indicates a partial derivative with respect to the coordinate that follows. 

Generally, when a beam is subjected to arbitrary loading, the normal stress in its longitudinal (axial) direction is 
much larger than those in its transverse directions. The stress state at an arbitrary point in the beam should be as 
shown in Fig. 1 (b), a space stress state with zero normal stress cr yy and o zz as well as zero shear stress r . 

Fortunately, the strain state is matched, as shown in equation (5), s yy = s zz = y yz — y zy — 0 . 


Also without loss of any generality, let us consider that the beam shown in Fig. 1 (a) is subjected to two 
concentrated loads, P y and P z , in the y and z directions, respectively, and a concentrated couple, P e , around the 

axis x, at the top end of the beam. This problem becomes a spatial cantilever beam subjected to static loads. Using 
equation (5) and Hook's law between stresses and strains of an elastic member with small deformations, the elastic 
potential energy of the beam is given by 


L 1 

n = J dx fj - (<T„e„ +T xy r xy + )dydz , 

0 O Z 


(6) 


in which, Q is the cross-sectional area in the yoz plane. 

The work done by non-conservative forces is 

Kc = PyVoy (L) + P z v 0z (L) + P g 0{L) , (7) 

where L is the length of the beam. 

The governing equations under these loading conditions will be obtained by resorting to the minimum potential 
energy condition, i.e. 8 (n - W nc ) =0 . This leads to the following. 

A W "(x) -BW(x)-CV \x) = 0 
DV"(x) + C t W'(x)= 0 



In which, 0 is a zero array, and 
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The boundary conditions corresponding to equation (8) are 

(W (0) = 0 
at x = 0, < , 

[V ( 0 ) = 0 


(9) 
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f AW'(L) = 

x = L, T . (10) 

[DV'(L) + C t W(L)-P = 0 


Equation (8), as well as its corresponding boundary conditions, (9) and (10), constitute a boundary value problem 
of ordinary differential equations, and its numerical solutions can be obtained by an Ordinary Differential 
Equation Solver (ODE Solver) [6-7] , which is of high quality and high efficiency, known as the numerical 
analytical solution. 


Numerical Examples and Results 



>K > | 


0.6m 0.3m 

FIG. 2 AN L-SHAPED BEAM AND ITS CROSS SECTION 


Figure 2 shows the cross section of an L-shaped cantilever beam of a length l = 3.0m . The Young's modulus of the 
material E = 3.0xl0 10 Pa and the shear modulus G = 1.22xl0 10 Pa . Other dimensions of the cross section are shown 
in the figure. In order to investigate the mechanical caracteristics of non uniforme torsion, non-symmetric bending, 
and coupling action of non-symmetric bending and torsion for the special thick-walled open section beam, 
respectively, three loading conditions are considered. 

Loading Condition 1 : a concentrated torque P e =500kN-m is applied at the free end; Loading Condition 2 : a 
transverse concentrated force of lOOkN towards the z direction is applied and passed through its geometric center 
at the free end; Loading Condition 3 : a force of 200kN in the y direction and a force of lOOkN in the z direction 
passed through its geometric center, and a counter clockwise torque of 500kN • m around the x axis are loaded at 
the free end. 

By means of the FNLM, the governing equations and corresponding boundary conditions of the L-section beam are 
identical to equtions (8) to (10) for the three loading conditions, except that the loading array in eqution (10) will 
respectively become 

Condition 1: P = [0.0 0.0 500.0xl0 3 f, 

Condition 2: P = [0.0 100.0xl0 3 O.Of, 

Condition 3:P = [200.0xl0 3 100.0xl0 3 500.0xl0 3 f. 

The boundary-value problem formulated by equations (8) to (10) with the above loading condition can be solved 
by ODE Solver. The stiffness matrices A, B, C and D can be determined by the size and the elastic modulus of the 
cross section of the example; the centroid coordinates can be determined to be at position (0.3m from long side, 
0.45m from short side), needed for the determination of the stiffness matrices. 

Table 1 presents the normal stress at different position of nodal line A- A and warping displacements of different 
cross sections due to the restraint torsion; Table 2 shows deflections in the y and z directions, and torsional angles 
around the x axis (counter clockwise positive) at different cross sections of the beam due to the restrained torsion; 
Table 3 indicates deflections in the y and z directions, and torsional angles at different cross sections of the beam 
due to the transverse force towards the z direction. 
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Figures 3, 4, and 5 display the axial deformations of the cross-section at x=1.5m due to the non-uniform torsion, the 
non-symmetric bending and the coupling action of non-symmetric bending and torsion, which are respectively 
depicted by the FNLM numerical solutions. 

TABLE 1 THE WARPING STRESS AND WARPING DISPLACEMENT OF THE NODAL LINE A-A 
AT DIFFERENT SECTIONS DUE TO NON-UNIFORM TORSION 


x/m 

a (x)/MPa 

m(x)/ m 

0.00 

10.06 

0.00 

0.75 

1.92 

1.30xl0- 4 

1.50 

0.33 

1.52x10-4 

3.00 

0.00 

1.57x10-4 


TABLE 2 THE DEFLECTIONS AND TORSIONAL ANGLES OF L-SHAPED TABLE 3 THE DEFLECTIONS AND TORSIONAL ANGLES OF L-SHAPED 

BEAM AT DIFFERENT SECTIONS DUE TO NON-UNIFORM TORSION BEAM AT DIFFERENT SECTIONS DUE TO THE FORCE ALONG THE Z AXIS 


x/m 

vo,(*)/ m 

v°:(*)/m 

ff(x) /rad 

x/m 

V0y(x) /m 

vod*)/ m 

9(x) /rad 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.75 

-2.39x10-4 

1.46x10-4 

1.19xl0- 3 

0.75 

3.14x10-5 

7.18x10-5 

2.93x10-5 

1.50 

-6.61x10-4 

4.01x10-4 

3.01x10-3 

1.50 

1.19x10-4 

2.32x10-4 

8.02x10-5 

3.00 

-1.58xl0- 3 

9.51x10-4 

6.88x10-3 

3.00 

3.97x10-4 

6.93x10-4 

1.90x10-4 


torsional angle : 3.01 x 10 3 



FIG. 3 THE CROSS-SECTIONAL AXIAL WARPING 
DUE TO NON-UNIFORM TORSION 



FIG. 4 THE CROSS-SECTIONAL AXIAL WARPING 
DUE TO NON-SYMMETRIC BENDING 


unit: m torsional angle : 2.82x10 3 



FIG. 5 CROSS-SECTIONAL AXIAL WARPING DUE 
TO NON-SYMMETRIC BENDING AND TORSION 


Conclusions 

The maximum warping normal stress caused by pure restrained torsion will occur at the fixed end of the cantilever 
beam, while the maximum warping displacement induced by pure restrained torsion will appear on the cross 
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section at the free end of the beam, as shown in Table 1. 

When an L-shaped cantilever beam is subjected to restrained torsion, not only the torsional angles around x axis 
will be produced, but also the deflections both along the y direction and z direction will be generated as shown in 
Table 2. This indicates that an L-shaped thick-walled beam also has a shear center (torsional center), i.e., when a 
torque is not applied around the shear-centers connected line, both the torsional and bending deformations will be 
induced by the torque. 

When an L-shaped cantilever beam is subjected to a transverse force along the z axis, not only the deflection in the 
z direction will be yielded, but also the deflection in the y direction and torsional angle around x axis will be 
developed as shown in Table 3. This phenomena again shows that an L-shaped thick- walled beam has the shear 
center (or torsional center); the phenomena also demonstrates the characteristics of non-symmetric bending, i.e., 
when the load surface is inconsistent with the main centroidal inertial plane of the cross-section of a beam, the 
bending deflections will be produced in two directions. 

The cross section of a thick-walled beam with an aspect ratio less than 3 will be warping to a curved surface as 
shown in Fig. 3, 4 and 5, whenever it is under the action of a combination of bending and restrained torsion, 
restrained torsion alone, or non-symmetric bending alone. Especially, under restrained torsion alone, the cross 
sectional warping of an L-shaped cantilever is more obvious, and the warping normal stress due to restrained 
torsion (see Table 1) must be taken into consideration. 
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